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OVERVIEW 

In  this  report  we  describe  the  research  performed  at  The  Ohio  State  University 
on  the  project  Controller  Design  for  Unstable  Aeroelastic  Systems^  sponsored  by  the 
Air  Force  Office  of  Scientific  Research  (AFOSR),  between  June  15,  1993  and  January 
14,  1995,  under  contract  F49620-93-1-0288. 

The  main  objective  of  this  research  was  the  robust  controller  design  for  active 
flutter  suppression,  in  the  presence  of  time  delays  and  neglected  dynamics  in  the  feed¬ 
back  loop.  For  this  purpose  we  have  developed  an  algorithm  (and  implemented  on 
the  computer)  which  derives  a  set  of  finitely  many  linear  equations  and  computes  the 
optimal  and  suboptimal  controllers  for  unstable  systems  with  time  delays.  With 
this  algorithm  we  can  now  compute  the  effect  of  time  delays  on  the  best  achievable 
robustness  and  performance  level  in  unstable  aeroelastic  systems.  Moreover,  all  sub- 
optimal  controllers  are  computed  for  a  given  system  with  time  delay  and  dynamic 
uncertainty  profile.  Our  algorithm  determines  the  “rational  parts”  of  suboptimal 
controllers.  Such  controllers  contain  an  infinite  dimensional  fixed  part,  which  is  the 
given  time  delay.  Therefore,  the  controller  is  “infinite  dimensional.” 

The  problem  of  finding  a  finite  dimensional  suboptimal  controller  is  also  studied 
in  this  project.  For  the  SISO  case,  it  is  shown  that  by  approximating  the  infinite 
dimensional  parts  of  the  optimal  controller  (with  certain  interpolation  constraints),  a 
rational  suboptimal  controller  can  be  obtained. 

In  this  project  we  have  also  performed  robustness  analyses  of  controllers  designed 
for  active  flutter  suppression.  Several  different  robust  control  algorithms  are  compared 
in  terms  of  the  robustness  of  stability  with  respect  to  perturbations  in  the  air  speed 
and  unmodeled  small  time  delays.  Our  numerical  studies  on  a  thin  airfoil  example 
showed  that  the  best  results  are  obtained  when  weighted  gap  optimization  using 
optimal  control,  and  ^-synthesis  methods  are  employed.  Also,  studied  in  the  project  is 
the  computational  complexity  of  certain  robust  control  algorithms,  such  as  //-analysis 
and  synthesis,  simultaneous  stabilization,  and  solution  to  bilinear  matrix  inequalities 
which  appear  in  and  mixed  types  of  robust  control  problems.  We  have  shown 

that  many  of  these  problems  are  NP-hard,  and  hence  it  is  not  feasible  to  develop  a 
polynomial  time  algorithm  to  find  exact  solutions  to  these  problems. 

The  principal  investigator  of  this  project  was  Hitay  Ozbay.  One  doctoral  student, 
Onur  Toker  (expected  graduation  is  Summer  1995),  was  supported  by  this  grant,  and 
his  contributions  were  as  much  as  the  Pi’s. 

Four  journal  papers,  based  on  this  research,  are  accepted  for  publication  and  three 
others  are  under  review.  Three  of  the  accepted  papers  are  co-authored  by  H.  Ozbay 
and  O.  Toker,  and  one  of  them  is  single  authored  by  O.  Toker.  Preliminary  versions  of 
these  papers  already  appeared  or  accepted  for  publication  in  conference  proceedings. 
A  list  of  publications  resulting  from  this  research  are  given  in  the  last  section  of  this 
report. 
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1  Introduction 


The  main  goal  of  this  project  was  to  investigate  robust  control  algorithms  for  unstable 
aeroelastic  systems,  against  unmodeled  dynamics  (mainly  due  to  poor  modeling  in 
aeroelastic  loads)  and  in  the  presence  of  time  delays  in  the  feedback  loop. 

A  robust  control  algorithm  has  been  developed  for  general  MIMO  unstable  sys¬ 
tems  with  time  delays.  This  methods  follows  the  gap  metric  approach  for  uncertain 
systems.  Also,  the  general  theory  for  control  of  infinite  dimensional  SISO  plants 
have  been  simplified.  The  controller  expression  obtained  in  this  project  is  so  simple 
that,  computation  of  the  optimal  and  all  suboptimal  controllers  can  be  done  by  using 
finite  dimensional  linear  algebra,  see  [23],  [10].  Implementation  issues  for  such  con¬ 
trollers  have  been  discussed,  and  an  approximation  procedure  has  been  proposed  in 
[30],  |29]. 

Several  robust  control  algorithms  have  been  tested  for  active  flutter  suppression 
of  a  thin  airfoil,  see  [26].  In  this  study  it  has  been  shown  that  even  small  time  delays 
(not  considered  in  controller  design)  may  cause  instability,  or  poor  performance.  The 
MIMO  gap  metric  design  algorithm  developed  in  this  project  has  also  been  applied 
to  this  system.  The  results  show  closed  loop  stability  with  guaranteed  robustness 
measure,  with  respect  to  gap  metric.  But  the  controller  itself  contains  time  delays  and 
its  finite  dimensional  parts  are  very  large  order.  Nevertheless  the  Bode  plots  of  the 
controller  indicate  that  it  is  easy  to  find  a  low  order  finite  dimensional  controller  which 
will  do  the  job,  see  Section  4  below. 

Computational  issues  of  several  robust  control  algorithms  are  also  investigated; 
and  it  is  shown  that  certain  analysis  and  synthesis  problems  computationalh^  are  not 
feasible,  see  [27],  [25],  [24],  and  Section  5  below  for  precise  statements  of  the  results. 

The  rest  of  this  report  is  organized  as  follows.  In  Section  2  robust  control  prob¬ 
lems  for  aeroelastic  systems  are  reviewed.  In  Section  3  control  algorithm  for  MIMO 
time  delay  systems  is  presented.  Applications  of  different  robust  control  techniques  to 
active  flutter  suppression  are  given  in  Section  4.  In  Section  5  certain  results  concerning 
the  computational  issues  in  robust  control  are  summarized.  Concluding  remarks  are 
made  in  Section  6.  A  list  of  publications  based  on  this  project  can  be  found  in  the  last 
section. 
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2  Robust  Control  of  Aeroelastic  Systems 


The  most  common  approach  for  modeling  large  scale  structures  is  to  use  linear  finite 
element  (FE)  methods.  In  aeroelastic  systems,  where  there  is  coupling  between  the 
structure  and  the  air  flow,  dynamic  aeroelasticity  equation  can  be  written  as  (see  [1], 
[4],  [12],  [1.3],  etc.), 

Mx  Cx  +  Kx  =  F^,ero{x,x)  +  Fext  (1) 

where  a.’  is  the  FE  displacement  vector,  M,  C,  K  are  the  mass,  structural  damping  and 
stiffness  matrices  respectively;  Faero  and  Fext  represent  the  aerodynamic  forces  and  the 
external  (e.g.  control,  disturbance,  noise,  etc.)  forces  respectiveley. 

In  flutter  analysis  one  looks  for  singularities  of  the  equation  (1),  assuming  Fext  = 
0.  In  the  flutter  suppression  problem  the  purpose  is  to  choose  (if  possible)  Fext  as  a 
feedback,  i.e.  Fext(.^■)  is  an  operator  acting  on  x,  such  that  the  closed  loop  system 
is  stable.  The  key  difficulty  in  these  problems  is  the  computation  of  Faero,  which, 
in  theory,  should  be  an  infinite  dimensional  nonlinear  operator.  In  practice,  Faero  is 
expressed  as  an  output  of  a  linear  system  whose  input  is  x\  in  the  frequency  domain 
we  have 


Faero(jkt.’)  —  H (^J Lo) (^j Lo') 


(2) 


where  is  determined  from  the  geometry  of  the  structure,  and  the  aerod3aiamic 

conditions  (in  particular  the  air  speed  F),  see  e.g.  [5]  [6]  [19]  [12]  [1.3],  etc.  In  order  to 
determine  the  flutter  speed  we  need  to  analyze  the  location  of  the  eigenvalues  of  the 
matrix 

A(a)  :=  Ms'^  +  Ca  +  K  -  H{s) 

as  V  changes.  Flutter  occurs  at  the  air  speed  V  =  F/,  which  makes  A{ju!)  singular. 
In  other  words  Vj  is  the  critical  air  speed  at  which  the  aeroelastic  system  becomes 
dynamically  unstable. 

Since  H{s)  is  impossible  to  determine  precisely,  in  practice  an  approximate,  say 
//,i(s),  is  used  in  the  flutter  analysis,  and  suppression  problems.  Hence,  neglecting 
approximation  error  and  parametric  uncertainties,  which  may  occur  in  modeling  M, 
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Feedback  Controller 


Figure  1:  Aeroelastic  system  model  with  time  delay 

C,  K ,  and  assuming  zero  initial  conditions  and  taking  the  Laplace  transforms  of  both 
sides  of  (1)  one  gets 

X{s)  =  [ms^  -^08  + K-  (3) 

as  the  nominal  model  to  be  used  in  control  design,  and  stability  analysis.  For  feedback 
controller  design  purposes  we  can  think  that  Fe,^i{t)  =  5ou(f),  where  the  vector  u{t) 
represent  the  command  inputs,  which  are  the  outputs  of  a  controller,  whose  inputs  are 
the  measured  outputs.  The  mesured  outputs  are  represented  in  one  vector  y{t),  given 
by 


y{t)  =  C(ix{t)  +  u,>i(t). 


where  Wi{i)  is  the  measurement  noise.  The  noise  is  assumed  to  be  the  output  of  a  filter 
whose  input  is  a  finite  energy  signal  w{t).  Here  Bq,  jBj,  Cq,  are  appropriate  size 
matrices.  A  block  diagram  of  this  aeroelastic  system  model  is  shown  in  Figure  1,  where 
D{s)  represents  possible  time  delays  in  the  feedback  loop,  and  BiWn{s)  can  be  seen  as 
the  actuator  disturbance  or  an  artificial  noise  which  accounts  for  unmodeled  dynamics. 

In  the  literature,  several  different  control  schemes  for  active  flutter  suppression 
have  been  reported,  see  e.g.  [9]  and  references  therein.  But  none  of  these  controllers 
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Figure  2:  Standard  feedback  configuration 

take  into  account  the  unmodeled  dynamics  or  possible  time  delays  in  the  feedback  loop. 
In  the  following  section  feedback  controller  design  for  this  type  of  input  /output  models 
of  the  aeroelastic  system  will  be  discussed. 


3  A  Robust  Control  Algorithm  for  MIMO  Time 
Delay  Systems 


Now,  the  aeroelastic  plant  model  developed  in  the  previous  section  can  be  written  as 


P[s)  = 

where  =  D(s)  is  the  delay,  and 

a(.s)  =  CoiMs^  +  Cs  +  K  -  Hn{s)r^B„ 

is  a  strictly  proper  rational  m  x  n  transfer  matrix.  The  feedback  control  system  is 
shown  in  Figure  3,  (here  C{s)  is  the  controller,  and  P{s)  is  the  plant,  r  is  the  reference 
input,  d  is  the  disturbance).  The  closed  loop  system  [P,  C]  is  said  to  be  stable  if  the 
entries  of  all  transfer  function  matrices  (/  —  CP)”^,  (/  —  CP)~^C\  P{1  —  CP)“\  and 
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P{I  —  CP)  belong  to  H°°.  When  the  closed  loop  S3^stem  is  stable  we  can  define 


bp.c  = 


{I-CP)-^[I  -C] 


as  the  stabilitj-  robustness  level  of  the  system  [P,C].  Larger  the  bp^c  more  robustness 
we  have  against  unmodeled  normalized  coprime  factor  perturbations  of  P,  see  [8].  With 
this  definition,  closed  loop  systems  [Ps^C]  are  stable  for  all  Ps,  which  belongs  to  a  gap 
ball  of  radius  6  around  P,  if  and  only  if  ^  <  bp^c-  Hence,  for  a  given  plant  P,  the 
optimal  robustness  radius  can  be  defined  as 


bopi{P)  :=  sup  bpfi. 

C  stabilizes  P 

In  [18],  it  is  shown  that,  bopt{P)  can  be  computed  using  state  space  techniques.  The 
above  problem  is  called  the  gap  metric  optimization  problem,  and  has  another  interpre¬ 
tation  in  terms  of  the  (unweighted)  coprime  factor  perturbations  of  the  original  plant. 
Namely,  if  P  =  N is  a  normalized  coprime  factorization,  the  above  optimization 
problem  corresponds  to  finding  the  maximum  t  such  that  there  exists  a  controller  C 
which  stabilize  all  Pa  =  [N  -f  An){M  +  A^)"^  with  ||[  Ayv  Am]  ||oo  <  e.  This 
maximum  e  value  is  equal  to  bopt{P)- 

The  suboptimal  robustness  problem  is  to  find  a  parameterization  of  the  set 

C-y  =  |c  :  [P, C]  stable  ,  6p,c  >  7|.  (4) 

This  problem  is  equivalent  to  parameterizing  all  controllers  C  which  robustly  stabilizes 
P  against  all  coprime  factor  perturbations  of  the  form 

II  [  An  Am  ]  ||oo  <  7  =:  ^ 

VP^  +  1 

for  a  given  desired  robustness  level  7.  The  main  result  of  this  section  is  a  parameter¬ 
ization  of  obtained  in  terms  of  the  state  space  realizations  of  Po.  The  algorithm 
summarized  below  has  been  implemented  using  MATLAB  commands. 
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Outline  of  the  procedure: 


Let  {Ao,  B,  C)  be  a  minimal  realization  of  Po{s),  then  find  the  stabilizing  solution 
Bp  of  the  algebraic  Riccati  equation 

A,Rf  +  BpAl  -  RfC’CRf  +  BB*  =  0  (5) 

and  let  A  —  Ag  H C  where  H  =  —RpC*.  Then 

F{s)  ;=  [-Asm]  =  [-C{sl  -  A)-^Be-^%R^  +  C{sl  -  A)-'^H  . 

Similarly,  find  the  stabilizing  solution  Rq  of  the  algebraic  Riccati  equation 


A;Rg  +  RgAo  -  RgBB^Rg  +  =  0 


and  let  Aq  =  Ag-V  BHg  where  Hg  =  -B*Rg.  Then 


■In  +  Hoisl-Aar'^B] 
.  C{sl  -  Aoy^Be-'^^  . 


(6) 


Since  P  =  NM  ’  =  M  is  a  coprime  factorization  (in  fact  normalized  coprime 
factorization),  there  exist  U^V^tj^V  G  H°°  such  that  the  generalized  Bezout  equation 

i>  -in  iM  in  iM  u]  [  y  -U]  _  [/„  o  i 

-N  mJ[a^  v\  [iV  vj  [-A'  Af]~[o  /„J 

holds.  To  parameterize  the  set  of  all  suboptimal  controllers,  first  choose  a  positive  real 
number  a.  Then,  let  xl,yl  G  be  the  solution  of 


md  r:y:  =  - 


,  Im  1 

pXg  H - 1 - , 

p  s  +  a 


and  G  7^2  be  the  solution  of 


=  -pvo  and  T^yo  =  -px* 


In  1 

p  s  — 


(9) 
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where  the  Hankel  opei’ator  is  defined  as 


In  [22],  we  discussed  how  to  obtain  numerical  solutions  of  (8,9)  using  state  space 
techniques. 

Given  numerical  solutions  of  (8,9),  define 

=  [2aa-J(a)]-V2  »  q  =  [2a(ar2(a))  »  0 

PaaiA^^)  =  pA  +  PaakA)  =  pA  +  a)  {xl{s)^ 

Q\ai<A)  =  pA  ~  ‘^)(ylA)j  G\^j^-  QaakA)  =  pA  ~  <P)ylA)G\AK 

wdiere  >>  stands  for  positive  definite  square  root. 

Then,  the  set  of  all  suboptimal  controllers  can  be  parameterized  as  : 

C,  =  {Ki(s)  +  A^,,2(s)£:(s)][i2,,i(s)  +  Z),,2(s)f(s)]-i  :  SeH^  <  1), 

where 

KMA  =  I/-.W(5)Pi^„.(-s)  -  N-{s)Q\^gls)] 

DoAs)  =  l/)A'(.«)i3W(-«)  +  M-(*)rV(»)l 
DM  =  [pN^PhA-A -H(D(s)Q'‘AAKiA\- 

Note  that  all  suboptimal  controllers  are  expressed  as  linear  fractional  transformation 
of  £,  with  coefficients  determined  by  N,M,N,N'  and  solutions  of  (8,9). 

In  the  next  section  the  Bode  plots  of  a  controller  obtained  from  the  above  pro¬ 
cedure  is  presented  for  an  active  flutter  suppression  problem.  For  the  same  problem, 
results  of  different  robust  control  algorithms  are  also  given  in  the  next  section. 


9 


Two  dimensional  thin  airfoil 


Figure  3:  Two  dimensional  thin  airfoil,  ’o’  observation  points,  ’x’  actuator  locations 

4  Applications  to  Active  Flutter  Suppression 


The  purpose  of  this  section  is  to  give  a  comparison  of  different  robust  controller  design 
techniques  for  a  linear  model  of  a  two  dimensional  thin  airfoil  described  in  [7]. 


4.1  Mathematical  Model  of  the  Plant 

Consider  the  two  dimensional  thin  airfoil,  of  length  i  =  0.305m  and  with  c  =  0. 076777, 
shown  in  Figure  3.  In  this  system  two  point  force  actuators  are  located  at  x  =  3^/4, 
y  =  ±c/2,  and  vertical  displacements  w{x,y,t)  at  x  =  0.3^,  y  =  ±0.3c  are  assumed  to 
be  available  for  feedback  control. 

The  mathematical  model  of  this  system  is  taken  from  [7],  and  it  is  as  given  below. 


TrpcVC, 


Mi) 

dt 


+  [K  -  wpV^lQqit) 


-7:pcVD-ii]i{t)  -  TrpcVD27]2{i)  -  =  0,  (10) 
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dViji) 

dt 


2V 

-0.041 — hih{t) 
c 


A  <^^9(0 


H- 


■  c  ^  dt  ^ 


(h2{t) 

di 


-0.032— /5r72(0  - 

c 


4  <^^9(0 


+ 


y 

c  "  dt 


where  the  entries  of  2  x  1  vector 


(11) 

(12) 


F{t)  = 


[^1(01 

[^’2(0J 


represent  the  control  forces  (to  be  generated  bj^  the  feedback  controller),  the  entries  of 
•5x1  vector  q{t)  are  the  generalized  coordinates,  and  the  entries  of  5  x  1  vectors  T]\{t), 
i]2{t)  are  the  generalized  aerodynamic  coordinate  vectors.  The  nominal  value  of  the  air 
speed  (in  the  direction  of  j/-axis)  is  taken  to  be  V  =  Wmisec,  and  the  other  constants 
are 

p  =  1.225  kg/m® 
c  =  0.076  m 
£  =  0.305  m 


■  1.315 

1.315- 

0.270 

0.270 

II 

0.631 

-0.631 

-O.OSl 

-0.815 

.  0.800 

0.800. 

-  0.000 

0.000 

0.289 

-0.015 

-0.444- 

0.000 

0.000 

0.015 

0.294 

-0.541 

Ka  = 

0.000 

0.000 

0.069 

0.000 

0.000 

0.000 

0.000 

0.000 

0.072 

0.000 

.  0.000 

0.000 

0.037 

0.042 

0.681 . 

■-0..305 

0.000 

0.144 

-0.007 

-0.370- 

0.000 

-0.305 

0.007 

0.147 

-0.451 

-0.072 

-0.037 

0.000 

0.000 

0.186 

0.004 

-0.074 

0.000 

0.000 

0.209 

.-0.037 

-0.045 

-0.009 

-0.010 

-0.114. 
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£»i  =  -0.165/5,  D2  =  -0.355/5, 


■0.0283 

0.0000 

0.0000 

0.0000 

0.0000- 

0.0000 

0.0283 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0021 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0021 

0.0000 

.0.0000 

0.0000 

0.0000 

0.0000 

0.0034. 

•0.3050 

0.0000 

0.0000 

0.0000 

0.0190- 

0.0000  - 

-0..30.50 

0.0000 

0.0000 

0.0230 

Ma  = 

0.0000 

0.0000  - 

-0.0090 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000  - 

-0.0090 

0.0000 

- 

0.0190 

0.0230 

0.0000 

0.0000  - 

-0.1140. 

-  0.0121 

0.0000 

-0.0046 

-0.0056 

0.0328  - 

0.0000 

0.4749 

0.03.32 

-0.0573 

-0.2502 

K  =  lO'* 

-0.0046 

0.0332 

0.0326 

-0.0452 

-0.0691 

-0.0056 

-0.0573 

-0.0452 

0..3012 

0.0268 

.  0.0328 

-0.2502 

-0.0691 

0.0268 

3.3105  . 

A4  = 


4>i 

(i>2 

\4>^ 

6^5  J 


[<i>\  h  -\h  |<^5], 


^5  — 


<^1 

(f>2 

\4'3 

\<f>4 


6 


0  0  <j>3  (^4  —2(j>5 


=  cosh(^3;)  —  cos(^)  —  0.7340955 


sinh(yar)  -sin(— a-) 


(f>2  =  cosli(yx)  -  cos(y )  -  1.01846644 


.  1  /  “^'2  \  .  /  ^^'2  \ 

smh(y.T)  -sm(yaT) 
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(/)3  =  sin(5fi.r)  -  ^sinh(/i.T)  + 

Ji 


( qI  sin(^i^)  +  ^2/1  sinh(/if)\ 
V  9I  cos((;i^)  +  /i^  cosh(/2^')  / 


^cosh(/ia;) 


cos(^ia: 


(J2 


64  =  sin(c/2.'r)  -  —  sinh(/2a’)  + 

h 


i ff2  +  ^2/2  sillh(/2^)  \ 

V  92  cos(if20  +  fi  cosh(/20  / 


(^cosh(/2.T)  -  cosh(i^2a;)^ , 


?ii5  =  1  -  cos(— a;), 

where 

ni  =  1.8751041 
7?2  =  4.69409113 
/i  =  9.193529 
/2  =  17.022282 
5^1  =  6.373793 
g2  =  15.679999 

The  measured  outputs,  yi  and  1/2,  are  vertical  displacements  at  the  observation  points 
indicated  by  'o'  in  Figure  3,  i.e. 

yiit)  =  w{0.Si,-0.Zc,t),  y2{t)  =  zn(0.3£,  0.3c,  i) 


where 


w{a\y,t)  =  Ylyi{x,y)qi{t), 


7,(.T,?/)  =  (f)i{x)e,{y) 

^liy)  =  £2(2/)  =  1 

£3(y)  =  £4(2/)  =  y/c 

£5(y)  =  iy'^lc^  -  1/3 


In  general,  for  point  force  actuators  located  at  (a^i,  j/i ),...,  (cCp,  j/p),  the  matrix  F  is 
given  b}^ 


F  = 


7i(-iu,yi) 


L7i(a’f.>2/p) 


75(a;i,yi) 


75 Vp)  ^ 
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4.2  Flutter  Comparisons 


Now  consider  the  plant  described  above  with  transfer  matrix  Po{s),  which  depends 
on  V'  (in  this  section  we  assume  that  there  is  no  time  delay  is  the  system).  The 
open  loop  system  and  the  closed  loop  system  stability  is  analyzed,  and  the  results 
are  shown  in  the  following  figures.  We  see  that  for  the  open  loop  system  the  flutter 
speed  is  Vj  =  17. 2m /sec.  We  designed  controller  for  the  nominal  air  speed  V  = 
30m/sec,  which  means  that  the  open  loop  system  is  unstable  at  this  operating  point. 
The  figures  5-9  illustrate  stability  robustness  with  respect  to  perturbations  in  V,  for 
controllers  designed  using  LQG,  unweighted  and  weighted  gap,  mixed  sensitivity, 
and  //  s3'nthesis  techniques.  The  results  can  be  summarized  as  follows: 


controller 

admissible  V  (m/sec) 

0  (open  loop) 

0.00  <  F  <  17.20 

LQG 

5.37  <  y  <  30.49 

gap  (unweighted) 

29.51  <  F  <  30.22 

gap  (weighted) 

11.89  <  F  <  40.01 

(mixed  sensitivity) 

1 0.35  <  F  <  40.15 

fi  synthesis 

10.27  <  F  <  40.24 

The  optimal  controller  guarantees  robust  stability  for  all  plants  of  the  form 
P{s)  =  Po(s)  +  Aa(s),  where  the  uncertainty  is  bounded  by  ra{j(jo)  >  ||Aa(iu;)||, 
i.e.  ra{j(.o)  can  be  seen  as  a  measure  of  the  stability  margin.  For  the  flutter  suppres¬ 
sion  problem  considered  here.  The  plot  of  the  largest  allowable  additive  uncertainty 
magnitude,  raijco),  versus  uj  is  shown  in  Figure  10,  for  two  different  time  delays  1msec 
and  10msec,  for  three  different  nominal  V,  ZOm/sec,  50m/sec  and  200m/sec.  As  the 
time  delay  and  V  increases  the  robustness  level  decreases  almost  uniformly  in  the  fre¬ 
quency  spectrum.  The  additive  uncertainty  is  allowed  to  be  large  in  the  frequency 
range  Ird/sec  to  lOOrd/sec. 


4.3  Effects  of  Time  Delays 


Consider  the  controllers  Ch“  and  Clqg  obtained  from  the  and  LQG  designs  when 
a  time  delay  of  1msec  is  neglected.  Time  domain  responses  (the  initial,  t  =  0,  and 
the  final,  t  =  4.0sec,  shapes  of  the  structure)  to  a  unit  gust  disturbance  are  shown  in 
Figures  11  and  12,  for  the  closed  loop  system  with  no  delay  and  1msec  delay.  We  see 
that  the  (LQG)  controller  is  (is  not)  robust  to  unmodeled  time  delay  of  1msec. 
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velocity  (m/sec) 


Figure  4:  Flutter  analysis  (open  loop) 


Flutter  Speed  analysis  (LOG) 


Figure  5:  Flutter  analysis  (LQG) 


Flutter  speed  analysis  (GAP) 


Figure  6:  Flutter  analysis  (unweighted  gap) 


Flutter  speed  analysis  (WGAP) 


Figure  7;  Flutter  analysis  (weighted  gap) 
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velocity  {m/sec) 

Figure  8:  Flutter  analysis  mixed  sensitivity  minimization) 


Flutter  speed  analysis 


Frequency  (rad/sec) 


Figure  10:  Largest  Allowable  Additive  Uncertainty 

Of  course,  it  is  better  to  include  the  time  delay  in  the  plant  model,  if  it  is  known. 
In  this  case  a  controller  can  be  designed  to  achieve  better  performance  and  stability 
properties.  Optimal  robustness  level  bopt,  corresponding  to  this  aeroelastic  system  (with 
transfer  matrix  P{$)  =  e“^U000Po('S),  where  Pp  depends  on  the  air  speed  V)  for  three 
different  values  of  V,  as  a  function  of  the  time  delay  T  is  shown  in  Figure  13.  Since  the 
unweighted  gap  does  not  give  a  reasonable  robustness  with  respect  to  perturbations  in 
V  we  chose  to  design  a  weighted  gap  for  the  original  plant  scaled  by  1000.  Note  that 
time  delaj's  in  the  order  of  10msec  leads  to  poor  robustness.  That  is,  bopt  is  considerably 
smaller  when  delay  is  around  10msec  or  higher.  Whereas  time  delays  up  to  0.2msec  do 
not  affect  the  optimal  robustness  level  signiflcantlj^  The  same  figure  also  shows  that 
as  V  increases  the  optimal  robustness  level  bopt  decreases. 

For  the  numerical  values  V  =  30?n/sec  T  =  1msec,  we  find  that  optimal  weighted 
gap  robustness  is  bmax  =  0.32.  If  we  choose  bsub  =  0.30  as  the  suboptimal  robustness 
level,  and  apply  the  controller  design  algorithm  we  developed  for  MIMO  unstable  delay 
sj'stems  we  obtain  a  controller  C{s)  which  gives  bp^c  =  0.30.  Note  that  in  this  numerical 
example  the  controller  is  a  2  x  2  transfer  matrix  of  the  form 

C(s)  =  {Dn  +  B„{sl-An)-'^Cn)({Dd  +  BdisI-Ad)-'^C,) 

+e  ^^{Ddd Bdd{sl  —  Add)~^Cdd)J 

The  MATLAB  based  program  written  for  this  problem  generates  the  matrices 
Bni  Ani  Cn,  Dd,  Bd,  Arf,  Crf,  Ddd^,  Bdd,  Addi  Cdd-  In  this  example  it  turns  out  that 
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t=0.0  sec,  Delay=0  msec 


t=4.0  sec,  Delay=0  msec 


t=0.0  sec,  Delay=1  msec 


t=4.0  sec,  Delay=1  msec 


Figure  11:  Initial  and  final  shape  of  the  structure  with  H°°  control 


t=0.0  sec,  Delay=0  msec 


t=4.0  sec,  Delay=0  msec 


Figure  12:  Initial  and  final  shape  of  the  structure  with  LOG  control 


the  matrices  A,,,  Ad  and  Add  have  dimensions  294  x  294,  236  x  236,  and  114  x  114 
respectivel}^  Therefore,  the  finite  dimensional  parts  of  the  controller  is  very  high  order. 
However,  the  magnitude  and  phase  plots  of  the  entries  of  the  controller  suggest’  that 
there  should  be  a  low  order  approximation.  At  this  point  any  approximation  algorithm 
can  be  used  (see  [11]  and  references  therein)  to  implement  a  low  order  controller. 


5  Remarks  on  the  Computational  Aspects 


This  section  summarizes  the  results  obtained,  during  the  project,  on  computational 
aspects  of  certain  robust  control  problems. 


5.1  control  of  SISO  infinite  dimensional  plants 


Consider  the  closed  loop  system  shown  in  Figure  15,  where  P{s)  and  C{s)  are  the 
transfer  functions  of  the  plant  and  the  controller  respectively.  In  this  section  we  give 
a  very  simple  expression  for  the  optimal  and  suboptimal  controllers.  These  con¬ 
trollers  are  computed  by  solving  a  finitely  many  linear  equations,  which  can  be  written 
directly  form  the  controller  structure. 

It  is  assumed  that,  the  plant  has  finitely  many  unstable  poles  denoted  by  Qi,  . . . ,  at, 
and  the  transfer  function  can  be  factored  as 


mn{s]No{s) 

m.dis) 


where 


/ 

md{s) = n 


/i:=l 


^  -  Q  a- 

S  +  Ok  ’ 


lUn  €  H'^  is  inner  (i.e.  all-pass  function)  possibly  infinite  dimensional,  and  N^,  € 
is  outer  (i.e.  minimum  phase)  possibly  infinite  dimensional.  It  is  also  as.sumed  that 
0!i,. .  .,0!i  e  C+  are  distinct.  Let  5  =  (1  +  PC)~'^  and  T  =  1  -  5,  and  Wi,  W2  be  given 
weighting  functions.  The  optimal  mixed  sensitivity  problem  is  to  find 


To 


inf 

C  stabilizes  P 


(13) 


and  the  optimal  controller,  denoted  by  Copt^  corresponding  to  the  optimal  perfor- 


* 


22 


20 


Magnitude  C_1 1 


V 


Figure  15:  SISO  Feedback  System 


mance  level  7^.  That  is,  Copt  stabilizes  the  plant  P,  and  yields 

r  Wj(i  +  pCopty^  1 

[W2PCopf{l  +  PCopt)-^\ 

The  suboptimal  mixed  sensitivity  problem  is  to  parameterize  the  set 


Cp  = 


C  stabilizes  P, 


H/jP 

W2T 


00 


<P 


(14) 


for  a  given  desired  performance  level  p  >  Jq. 

In  this  section,  it  is  assumed  that  1^1(5)  and  H4(s)  are  rational  functions.  For 
properness  of  the  optimal  controller  it  will  be  assumed  that  ITifs)  is  non-constant  and 
ITi,  (W2K)  and  (fF2A^o)-'  e 

Let  rii,. . .  ,7]^^  e  C+,  be  the  poles  of  lTi(— 5)  (if  rji  has  multiplicity  then  it  is 
assumed  to  be  repeated  ii  times  in  this  list),  and  set 


Ey{s)  := 


/TU(-.s)II/i(s) 

V  7^ 


(15) 


The  zeios  of  are  denoted  by  /9i, . . . ,  /?2„, .  Suppose  they  are  distinct  for  a  chosen  7. 
Then,  /3ps  can  be  enumerated  in  such  a  way  that  A, . . . ,  0^,  are  in  C+,  and  B,,.,  =  -3, 
JN'ow  dehne 


E^{s)  :=  G^{s)  J] 


A:=l 


rtk  -  'S 
rjk  +  5 


(16) 
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where  is  minimum  phase  and  determined  from  the  spectral  factorization 


G'-,(s)G'-y(— s)  :: 


1  -  ( - - 1)  E^{s) 

\  I  } 


(17) 


Then  (under  certain  genericity  assumptions-see  [28])  the  optimal  controller  is  given 
by 


Copi[»)  =  Ey^{$)m.d{s) 


N.(s)-^F,ML{s) 

1  +  mn{s)E^^{s)L{s) 


(18) 


where  L{s)  —  L2is)/Li{s),  and  Li{s),L2{s)  are  polynomials  of  degrees  less  than  or 


equal  to  (7?i  +  /  —  1)  and  they  satisfy 

0  =  IilA-)  +  m.n{/Sk)F-,.^{l3k)L2{0k)  k=l,...,ni  (19) 

0  =  Ti(Q:fc)  +  mniak)Fy^{ak)L2{ak)  k  =  (20) 

0  =  L2i-l3k)  +  mn{fik)F^^{/3k)Li{-/3k)  k=^l,...,ni  (21) 

0  =  L2{-ak)  +  mn{ak)F~,.^{Qk)Li{-ak)  k  =  l,...J.  (22) 


In  general,  if  7^  is  replaced  by  a  variable,  say  7,  in  equations  (19-22),  then  a  new  set 
of  linear  homogeneous  equations  is  obtained,  in  terms  of  2(721  +  0  unknown  coefficients 
of  Li  and  L2,  for  each  fixed  7.  In  [28]  it  was  shown  that,  under  certain  genericity 
assumptions,  7^  is  the  largest  value  of  7  for  which  there  is  a  non-trivial  solution  to 
these  2(?2i  +  /)  linear  homogeneous  equations.  That  is,  7^  can  be  found  by  plotting 
smallest  singular  values  of  the  matrix  representation  of  these  equations,  as  7  varies  in 
an  interval.  The  largest  value  of  7  for  which  the  plot  shows  a  zero  is  70.  It  w’as  also 
shown  that,  see  [28],  in  the  optimal  case  L2{s)  =  Ti(— s)  and  hence  the  number  of 
unknown  coefficients  in  these  coefficients  reduce  to  ?2i  +  /;  and  these  can  be  determined 
from  (19,20). 

Similarly,  all  suboptimal  controllers  are  in  the  form 


Csubopt{s)  =  Ep{s)mi{s) 


K{s)-'^Fp{s)Lu{s) 

1  +  'l'>'>'n{s)Fp{s)Lu{s) 


(23) 


where 


L2{s)  +  Lri-s)U{s) 
Li{s)  +  L2{-s)U{s) 


U  eB, 
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and  Li{s),L2{s)  are  polynomials  of  degree  <  ni  +  /  satisf5dng  (19-22)  with  7^  replaced 
by  p,  and  the  following  two  conditions: 

0  =  l2(-fl)  +  {Ep{a)  +  l)Fp{a)mn{a)Li{-a)  (24) 

1  =  Lii-a),  (25) 

for  some  arbitrary  a  E  R,  a  >  0,  a  ^  Oi’s  and  a  ^  for  all  i  =  1,...,/  and 
j  =  1, . . .  77i.  In  (24,25)  G  is  a  free  parameter.  For  different  values  of  a  one  can  obtain 
different  parameterizations  of  the  suboptimal  controllers. 


5.2  Finite  Dimensional  Controllers  via  Approximations 


The  structure  of  H°°  controllers  (optimal  and  suboptimal)  is  given  in  the  previous 
section,  and  it  is  shown  that  for  distributed  systems,  the  optimal  controller  is 
not  rational.  Therefore,  the  question  of  obtaining  a  stabilizing  finite  dimensional  H°° 
controller,  whose  performance  is  “close”  to  optimum,  arises  naturally.  If  a  controller 
satisfies  this  requirement,  then  it  is  said  to  be  “approximately  optimal.”  One  can 
approximate  the  infinite  dimensional  parts  of  the  optimal  controller  to  obtain  such  a 
controller,  and  this  method  has  been  studied  before  [17]  for  stable  plants.  For  un¬ 
stable  plants  it  was  shown  that  [14]  if  coprime  factors  of  the  optimal  controller  are 
approximated  separately,  then  under  mild  conditions,  the  resulting  finite  dimensional 
controller  is  approximately  optimal.  However,  computation  of  the  coprime  factors  of 
the  optimum  controller  requires  construction  of  a  complicated  system  of  linear  equa¬ 
tions,  (compare  the  formulae  of  the  previous  section,  with  the  ones  in  [16,  14]). 

One  possible  approach  to  obtain  a  finite  dimensional  controller  Cj  from  the  con¬ 
troller  C  given  in  (18)  is  to  approximate  the  infinite  dimensional  parts  of  the  controller, 
and  m„,  by  rational  functions,  N~j  and  7?7„/,  and  obtain 


E^^{s)md{s) 


N;/{s)F,Ss)Lis) 

1  +  mnf{s)F^^{s)L{s)‘ 


(26) 


Define 


Sf  =  {l  +  PCj)-\  Tf  =  l-  Sf,  Sopt  =  (1  +  PCoptr\  T,pt  =  1  - 


F-y.L 


1  4-  ninF^^L 


{ninf  -  ?77„),  An^  =  1  -  N^jNo. 
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The  following  lemma  gives  a  set  of  sufficient  conditions  for  the  closed  loop  system 
(P,  C/)  to  be  stable. 

Lemma  1:  The  closed  loop  system  (P,  C/)  is  stable  if 


(i)  m„f  interpolates  m„  at  a/s,  in  C+,  and  juj  axis  zeros  of  +  mnF^L2, 


(ii) 

sup 

R 

(iii) 

sup 

(lO'olT2"'(i^)ANp(icu)|  +  l7oTT2-'(icu)A^,7(ia;)A^,(ia;)A„,„(ia>)|)  <  1, 

u;6R 

(iv)  1  - 

n:}n 

0  € 

Now.  let  us  consider  the  performance  of  the  closed  loop  system  {P,Cf). 

Lemma  2:  The  following  conditions  guarantee  that  the  performance  7(P,  C/)  of  the 
clo,sed  loop  system  (P,  C/)  will  be  between  %  and  7^  +  e: 

(i)  sup  |W7'(,a,)A„,{j<u)|  <  ^(1  -  e), 

(ii)  sup  \Wf\joj)Am„{juj)\  <  -^(1  -  e), 

(iii)  sup  \W2\juj)AN,{juj)\  <  -^(1  -  e),  and 

ueR  472 

(iv)  sup \W^\jLo)Amn{j^)\  <  T^(l  -  e), 
u;GR 

where  e  is  a  positive  real  number  less  than  1. 

Note  that,  since  strictly  proper,  approximation  conditions  on  m.nj  can  be 

interpreted  as  a  uniform  approximation  in  the  low  frequency  range  subject  to  certain 
interpolation  conditions  in  the  closed  right  half  plane.  Interpolating  both  the  values  of 
???„  and  its  derivatives  at  ^iS  gives  better  7nn/’s  functions  which  is  likely  to  decrease 
the  value  of  the  expression 

{mnjjuj)  -  mnf{joj))F^{ju})L{juj) 
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around  ^{’s.  Note  that  the  supremum  of  this  expression  over  ju)  axis  appears  as  one  of 
the  conditions  of  Lemma  1. 

We  refer  to  [23]  and  [29]  for  a  numerical  example  illustrating  the  approximation 
procedure  proposed  in  this  section. 


5.3  On  the  NP-hardness  of  certain  control  problems 

In  this  section  we  simply  state  the  main  results  of  the  papers  [27]  and  [25],  on  the 
computational  complexity  of  certain  robust  control  problems.  The  reader  is  refered  to 
the  full  versions  of  these  papers  for  complete  details. 

In  [27],  it  is  shown  that  for  a  given  complex  matrix  M,  and  a  purely  complex 
uncertainty  structure  A,  the  problem  of  checking  whether  the  inecjuality 

(min{(T„i„,,.(A)  :  A  G  A,  det(/  -  M A)  =  0})"^  =;  fJA{M)  <  1 

holds,  is  j\l  P-hard.  It  is  also  shown  that,  the  problem  of  checking  whether  the  frequency 
domain  //,  ||M(s)||^,  of  an  LTI  system,  Af(s),  is  less  than  1,  and  the  problem  of  checking 
whether  the  best  achievable  fj,, 

of  a  linear  fractional  transformation  (LFT),  is  less  than  one,  are  both  MV- 

hard  problems.  In  other  words,  purely  complex  //  computation,  analysis/synthesis  are 
A/’T’-hard. 

It  w'as  known  that  the  computation  of  is  A^P-hard  for  purely  real  [20]  and  mixed 
real/complex  [3]  uncertainty  structures,  but  these  results  do  not  give  much  information 
about  the  computational  complexity  of  the  purely  complex  //  problem,  nor  they  imply 
much  about  complexity  of  the  n  analysis/synthesis  problems.  Although  general 
norm  computation,  analysis/synthesis  have  a  well  established  theory  for  LTI  systems, 
there  is  no  known  non-conservative  polynomial  time  procedure  for  purely  complex  p 
computation,  analysis/synthesis  problems.  In  this  part  of  the  project,  [27],  we  gave 
proofs  of  the  A/'T’-hardness  of  the  above  mentioned  problems.  These  results  imply 
that  it  is  rather  unlikely  to  find  non-conservative  polynomial  time  procedures  for  the 
purely  complex  p  computation,  analysis/synthesis  problem,  contrary  to  the  standard 
H°°  problems. 
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As  independent  results,  it  is  also  shown  that,  [27],  the  problem  of  checking  the 
stability  and  the  problem  of  computing  the  norm,  are  both  A/”P-hard  problems 
for  multidimensional  systems.  These  results  imply  that  it  is  rather  unlikely  to  find 
a  simple  analogue  of  the  Schur-Cohn  test  for  checking  the  stability  and  an  efficient 
generalization  of  bisection  method  of  [2]  for  computing  the  norm,  in  the  context 
of  multidimensional  systems. 

A  matrix  valued  function,  F{X,  Y),  is  called  bilinear  if  it  is  linear  with  respect  to 
each  of  its  arguments.  An  inequality  of  the  form  F{X,  Y)  >  0  is  called  a  bilinear  matrix 
inequality  (BMl).  Recently  it  has  been  shown  that  many  robust  control  problems  can  be 
transformed  to  finding  a  solution  to  bilinear  matrix  inequalities  (BMIs).  For  example 
the  static  output  feedback  problem  and  fixed  order  robust  controller  design  problem 
can  be  put  in  this  framwork.  In  [25]  we  have  shown  that  the  problem  of  checking  the 
solvability  of  a  given  BMI  is  NP-hard.  Although  BMI  approaches  to  robust  control 
seems  to  be  a  potentially  powerful  tool,  it  is  rather  unlikely  to  find  a  polynomial  time 
algorithm  for  solving  general  BMI  problems.  In  [25]  it  is  also  shown  that  simultaneous 
stabilization  with  static  output  feedback  is  an  NP-hard  problem,  namely  given  n  plants, 
the  problem  of  checking  the  existence  of  a  static  gain  matrix  K  which  stabilizes  all  of 
the  n  plants,  is  NP-hard. 


5.4  On  the  order  of  simultaneously  stabilizing  compensators 

In  this  project  the  problem  of  simultaneous  stabilization  using  a  dynamic  compensator 
has  also  been  studied.  The  results  will  appear  in  [24],  where  it  is  shown  that  there  is 
no  upper  bound  for  the  minimal  order  of  a  simultaneously  strongly  stabilizing  com¬ 
pensator,  in  terms  of  the  given  plant  orders.  A  similar  problem  was  also  considered  in 
[21],  where  it  was  shown  that  such  a  bound  does  not  exist  for  the  strong  stabilization 
problem  of  a  single  plant.  But  the  examples  given  in  [21]  were  forcing  an  approximate 
unstable  pole-zero  cancellation,  or  forcing  the  distance  between  two  distinct  unstable 
zeros  to  go  zero.  In  [24]  it  is  shown  that:  (i)  if  approximate  unstable  pole-zero  cancel¬ 
lation  does  not  occur,  and  the  distances  between  distinct  unstable  zeros  are  bounded 
below  by  a  positive  constant,  then  it  is  possible  to  find  an  upper  bound  for  the  min¬ 
imal  order  of  a  strongly  stabilizing  compensator;  (ii)  and  for  the  simultaneous  strong 
stabilization  problem  (even  for  the  two  plant  case),  such  a  bound  cannot  be  found. 
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6  Conclusions 


Ill  this  project  robust  control  techniques  are  investiga.ted  for  unstable  aeroelastic  sys¬ 
tems.  A  new  algorithm  is  developed  (and  its  MATLAB  based  program  is  written) 
for  MIMO  delay  systems  to  optimize  robustness  measured  in  the  gap.  The  algorithm 
can  easily  be  modified  to  obtain  controllers  from  weighted  gap  optimization  method. 
Numerical  simulations  have  shown  that  when  time  dela3^s  are  neglected  in  the  design 
of  LQG,  gap,  loop  shaping  via  weighted  gap,  H°°,  and  ^  controllers,  the  closed  loop 
stability  may  be  lost  or  the  performance  may  be  poor  as  far  as  flutter  suppression 
problem  for  a  thin  airfoil  is  concerned.  Our  algorithm  has  been  tested  on  the  same 
example,  and  an  approximable  infinite  dimensional  controller  is  obtained. 

In  the  SISO  case  numerical  computation  of  controllers  is  simplified  consid- 
erabh'.  A  simple  MATLAB  based  program  is  written  to  compute  the  optimal  perfor¬ 
mance,  optimal  controller,  and  all  suboptimal  controllers.  Also,  a  numerical  approxi¬ 
mation  procedure  is  developed  to  derive  finite  dimensional  suboptimal  controllers  for 
general  infinite  dimensional  plants. 

Computational  complexity  of  certain  control  algorithms  are  studied.  In  particular 
it  is  shown  that  purely  complex  fi  computation,  analysis  and  synthesis  problems  are 
NP-hard;  and  so  are  the  problems  of  simultaneous  stabilization  with  static  output 
feedback  and  solution  to  bilinear  matrix  inequalities.  It  is  important  to  note  that 
the  above  problems  have  been  considered  by  many  researchers  in  the  past,  and  no 
polynomial  time  algorithmic  procedure  has  been  developed  for  exact  computation.  The 
results  of  this  research  indicate  that  it  is  practically  impossible  to  find  a  poljmomial 
time  exact  computation  procedure  for  the  above  mentioned  problems.  So,  the  scientific 
community  should  be  satisfied  with  approximate  solutions,  which  are  computational!}' 
more  feasible. 

As  far  as  the  aeroelasticity  application  (flutter  suppression),  simulations  show 
that  based  controllers  give  good  robustness  and  performance  compared  to  more 
conventional  LQG  method.  Of  course  in  both  cases  weight  selection  is  a  major  issue. 
Even  though  we  now  have  tools  to  compute  H°°  and  LQG  controllers,  in  many  cases 
for  a  given  set  of  weights  and  the  plant  the  “central  controller”  is  unstable  (see  [15] 
for  a  aircraft  pitch  angle  command  tracking,  and  gust  alleviation  examples).  On  the 
other  hand  for  safety  reasons,  in  the  actual  implementation,  unstable  controllers  are 
undesirable.  Therefore,  it  would  be  a  major  contribution  if  one  develops  methods  to 
determine  stable  suboptimal  controllers.  The  PI  will  be  focusing  on  this  problem  in 
the  future. 
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7  List  of  Publications  Based  on  the  Project 


1.  0.  Toker  and  H.  Ozbay,  “Gap  metric  problem  for  MIMO  delay  systems:  Param¬ 
eterization  of  all  suboptimal  controllers,”  Automatical  to  appear. 

2.  0.  Toker  and  H.  Ozba}^,  Optimal  and  suboptimal  controllers  for  infinite 

dimensional  SISO  plants,”  IEEE  Transactions  on  Automatic  Control,  to  appear. 

3.  C.  Gu,  O.  Toker  and  H.  Ozbay,  “On  the  two  block  problem  for  a  class  of 
unstable  distributed  systems,”  Linear  Algebra  and  its  Applications,  to  appear. 

4.  0.  Toker,  “On  the  order  of  simultaneously  stabilizing  compensators,”  IEEE 
Transactions  on  Automatic  Control,  to  appear. 

5.  0.  Toker  and  H.  Ozbay,  “On  the  NP-hardness  of  solving  bilinear  matrix  inequal¬ 
ities  and  simultaneous  stabilization  with  static  output  feedbaek,”  Proceedings  of 
the  American  Control  Conference,  Seattle  WA,  1995,  to  appear. 

6.  0.  Toker  and  H.  Ozbay,  “Robustness  analysis  of  controllers  designed  for  active 
flutter  suppression,”  Proceedings  of  the  American  Control  Conference,  Seattle 
W.4,  1995,  to  appear. 

7.  0.  Toker  and  H.  Ozbay,  “On  the  NP-hardness  of  the  purely  complex  mu  compu¬ 
tation,  analysis/synthesis,  and  some  related  problems  in  multidimensional  sys¬ 
tems,”  Proceedings  of  the  American  Control  Conference,  Seattle  WA,  1995,  to 
appear. 

8.  0.  Toker  and  H.  Ozbay  “On  the  structure  of  controllers,”  Proceedings  of 
33rd  IEEE  Conference  on  Decision  and  Control,  Lake  Buena  Vista,  December 
1994,  pp.  3949-3954. 

9.  O.  Toker  and  H.  Ozbay,  “Numerical  issues  in  the  computation  of  II°°  con¬ 
trollers  for  distributed  plants,”  Proceedings  of  31st  Allerton  Conference  on  Corn- 
munication,  Control  and  Computing,  University  of  Illinois  Urbana-Champaign, 
September- October  1993,  pp.  535-544. 
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